Dimensional crossover and driven interfaces in disordered ferromagnets 
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We study the depinning transition of driven interfaces in thin ferromagnetic films driven by external 
magnetic fields. Approaching the transition point the correlation length increases with decreasing 
driving. If the correlation length becomes of the order of the film thickness a crossover to two 
dimensional behavior occurs. From the corresponding scaling analysis we determine the exponents 
characterizing the transition of the three dimensional system. 
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I. INTRODUCTION 

We study interface motion in the random-field Ising 
model. In this model an interface separates regions of op- 
posite spin orientation. A magnetic field forces the inter- 
face to move whereas this motion is hindered by the dis- 
order. At zero temperature a permanent interface motion 
is found if the driving field H exceeds a critical threshold 
He determined by the disorder. The vanishing of the in- 
terface velocity at {H = Hc\T = 0) is in general referred 
to as a continuous phase transition Considering 
three dimensional systems, for instance, the interface ve- 
locity V vanishes at the so-called pinning/depinning tran- 
sition according to v{T = 0) ~ {H — Hc)^^'' and v(H = 
He) ~ T'^/^3d^ respectively (see e.g. The correla- 

tion length diverges algebraically: ^ ^ {H~Hc)~''-^'^ iH'ii- 

This pinning/depinning transition can be understood 
within the framework of the rcnormalization group the- 
ory (see e.g. [|6[-p^ and references therein). Within this 
theory it is found that the order parameter of a contin- 
uous phase transition is a generalized homogenous func- 
tion of, in general, several thermodynamic parameters . 
In many cases temperature and magnetic field belong 
to these parameters. Additional parameters may cause 
crossovers. Examples are spatial anisotropics in Heisen- 
berg models dipolar effects or restrictions due 
to geometry like in thin films (see e.g. and references 
therein) . In the case of thin films the system behaves like 
a three dimensional system only as long as the correla- 
tion length ^ is small as compared to the film thickness 
I. Approaching the transition point the component of 
the correlation length perpendicular to the film layers is 
bounded by the film thickness. If becomes of the 
order of unity a crossover from three to two dimensional 
behavior occurs. Experiments which determine domain 
wall velocities in magnetic films typically image the mag- 
netic state of a sample by looking onto the sample using 
the magneto-optical polar Kerr effect and a CCD cam- 
era fl^Jl^ . If the film is sufficiently thin it is possible to 
obtain the interface position from snapshots generated by 
the CCD camera and to calculate the interface velocity 
from the time dependence of this position. This was done 
in where a Pt(3.4nm)/Co(0.5nm)/Pt(6.5nm) film 
with perpendicular anisotropy was investigated. Due to 
the film thickness the authors saw evidence to neglect the 
height dependence of the interface position and to treat 



the interface not as a two dimensional interface but as 
a one dimensional line. With increasing film thickness, 
however, this assumption fails because then the corre- 
lation length drops below the film thickness causing a 
crossover to three dimensional behavior. This scenario is 
investigated in the present paper in the context of driven 
interfaces in the random-field Ising model. 



II. SIMULATION 

The RFIM is defined by the Hamiltonian 

4 E'^'^j^^E^'-E'^^'^- w 



n 



The first term is the exchange interaction of neighbor- 
ing spins and the sum is taken over all pairs of them 
{Si — ±1). The second term specifies the coupling to the 
driving field H . Additionally, the spins are coupled to in- 
dependent quenched local random-fields hi characterized 
by their probability density p{hi) given by 



P{hi) = i 



for 



M < A 
otherwise. 



(2) 



We use a random-sequential update with transition prob- 
abilities according to a heat bath algorithm in the limit 
of zero temperature. Since in the vicinity of the criti- 
cal point finite-size effects may not only occur due to the 
finite film thickness but also due to the other finite exten- 
sions of our system, we calculated the interface velocities 
for each film thickness I varying the other extensions of 
the system. For sufficient large extensions we observed 
no finite-size effects from which we concluded that the 
data presented in the following correspond within negli- 
gible errors to those of an infinite extended film of thick- 
ness I. Note that for this analysis the extensions of the 
system must be increased if one approaches the critical 
point. But this increase requires an increase in CPU time 
which restricts therefore the data to values not too close 
to the critical point. We investigate system sizes of up 
to Z X = 8 X 1024^ unit cells of a body centered cubic 
lattice. 

We apply periodic boundary conditions in the direc- 
tions parallel to the film and antiperiodic ones perpen- 
dicular to the interface (see j^). The interface moves 
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FIG. 1. Snapshot of a moving interface obtained in a film 
consisting of 6 magnetic layers. Overhangs, which are rare 
and small, are not shown. In the gray sketched area be- 
low the interface the spins are aligned parallel to the driving 
field. Spins above the interface are aligned antiparallel and 
not shown. 



along the [100] direction of the bcc lattice resulting in 
a finite interface velocity for any driving field if 7^ in 
the absence of disorder. The same behavior is found, 
for instance, on simple cubic lattices with an interface 
moving along the diagonal direction of the lattice [^|J|] 
and on diamond lattices with the interface moving along 
the [100] direction. The schematic phase diagram in y] 
applies to all of these cases. In particular, a continuous 
phase transition is found for A > J as long as nucle- 
ation does not occur. This is the case for A = 3, which 
turns out to be a convincing choice because then the di- 
mensional crossover is numerically accessible for a broad 
range of system sizes. 



III. RESULTS 

In our simulations we start with an originally flat in- 
terface which is built into the system. After a transient 
regime the interface reaches a stationary state. In Fig. |l| 
we show the snapshot of an interface configuration (for 
_ff = 1.8 and A = 3.0). Overhangs are not displayed since 
they are rare and small. We obtain the interface veloc- 
ity from the mean interface position h at a given time t 
according to v — dh/dt. In [^ it was argued that this 
velocity is a generalized homogeneous function of tem- 
perature and driving field. We generalize this ansatz to 
the present situation and include additionally a thickness 
dependence of v, 



(3) 



Here, T denotes the temperature and -q = H — H^'^ 
the reduced driving field. In a bulk system we recover 
under the assumption v{l^^ = 0) = finite [Eq. (||)] 
y ^ T^/Ssd f (^ /X^/PsdSsd^i which is known to be sat- 
isfied in the RFIM [§,§. In the following we consider 
only T = 0. Choosing A = [-l^sd/i^sd Eq. (||) we obtain 
for the velocity the following scaling behavior: 



(4) 



In this limit the value of Hr 



with f{x 00) ^ a;' 
coincides with that of the bulk system. For any finite 
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FIG. 2. Interface velocities v obtained for various film 
thicknesses I and different driving fields H . We also plot the 
interface velocity of the bulk system. Approaching Hc{l) the 
system size must be increased in order to find negligible fi- 
nite-size effects (see Sec. ^). The curves terminate at fields 
where the finite size effects set in for the largest set sites ac- 
cessible numerically. 



however, becomes I dependent and is shifted accord- 
ing to Hc{l) - H^"^ « x^l-'^'"^'' (see ^ and references 
therein). Inserting this relation into Eq. (^ one finds 
that f{x) does not vanish at a; = but at a finite value 

X*. 

In Fig. H we plot v{H) for different film thicknesses. 
For a comparison, we also plot velocities obtained in a 
three dimensional system. From the data it is evident 
that the curves v{H) for different film thicknesses deviate 
more and more from the bulk behavior with decreasing 
H. In the corresponding region of if -values the crossover 
from 2 to 3-dimcnsional behavior occurs which thus is nu- 
merically accessible. In the crossover region the interface 
velocity in a film turns out to be smaller than in the bulk 
meaning that the threshold field He is shifted towards 
larger values of the driving field. Rescaling the veloci- 
ties according to Eq. (^) yields the data collapse shown 
in Fig. ||. From the collapse it may be concluded that 
the scaling function f{x) [see Eq. vanishes at a value 
X* > 0, again meaning that the threshold field is shifted 
towards larger values with decreasing film thickness. We 
obtain = 0.677 ±0.07 and iJ^d = 1.491 ±0.02. These 
values coincide within the error-bars with those found 
in bulk systems. In this case (inset of Fig. |3|) we find 
^buik ^ Q g4 ^ Q Q5 ^buik ^ 5 ^ Q confirming 

the results of the crossover scaling. From the data col- 
lapse in Fig. ^ we also obtain — 0.763±0.03. Both P^d 
and v^d coincide within the error-bars with values found 
on other lattices and different values of A Also, 
they agree with the exponents of the Edwards- Wilkinson 
equation with quenched disorder (see e.g. [^jlO 1^ and 
references therein). For this equation /3qew ~ 0.62 and 
'^QEW ~ 0.77 was obtained by an e-expansion valid to 
order [0 . 

Since the interface velocity satisfies the ansatz Eq. (0) 
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FIG. 3. Scaling plot according to Eq. 
collapse we obtain fi^d ~ 0.677 ± 0.07, 
and H^'^ = 1.491 ± 0.02. The inset shows the interface ve- 



From the data 
= 0.763 ± 0.03 



locities of the bulk system. We obtain (3 



0.64 ± 0.05 



and = 1.5 ± 0.015. For small values of 77 aU curves 

join the data collapse. For sufficient large 1% i.e. well above 
the crossover region, the scaling ansatz Eq. (W) does not hold 
as can be seen from the deviations of single curves from the 
scaling behavior. 



it is possible to draw conclusions about the critical be- 
havior at Hc{l). Considering f{x) with x ^ rj 

/iZ-^^d and 

taking into account that i>3d > one finds a; = at Hc{l) 
for any finite film thickness. On the other hand, a three 
dimensional system corresponds to the limit I — > 00 and 
in this case f{x) is determined by the limiting behavior 
f{x — ^ 00) x^^"^ . The critical behavior in a film is 
therefore different as compared to the bulk behavior and 
the critical exponents (3 and v coincide with those of the 
two dimensional model for any finite film thickness. In 
the two dimensional model [32d ~ 0.31 and V2d ~ 1 was 
found [§,1). Taking Eq. (||) into account, we also conclude 
that the exponent 6 characterizing the thermal rounding 
of the depinning transition [v ^ T^/^2d at H = Hc{l)] is 
given by S2d ~ 5 (see |^) for any finite film thickness. 



IV. CONCLUSION 

In conclusion, we have investigated the depinning tran- 
sition of a driven interface in thin films. We have found 
that the critical behavior is governed by the two di- 
mensional fix-point and the corresponding exponents for 
any finite film thickness. The exponents obtained by 
the crossover scaling are in agreement with those of 
the Edwards- Wilkinson universality class. The scaling 
ansatz used to analyse our data could also be used to 
determine critical exponents of the depinning transition 
of three dimensional systems experimentally, in particu- 
lar since present experimental techniques use thin films 
as samples |l^,|lj]. In jl3|, for instance, Co28Pt72 alloy 
films were investigated with grain sizes of typically 20 nm 
and film thicknesses of 5...50nm. If one naively assumes 



that by a variation of temperature and/or driving field it 
is possible to increase the correlation length from the size 
of the grain to the film thickness we expect the crossover 
scaling [Eq. (^)] to work and to yield the exponents of 
a three dimensional sample. Note that due to dipolar 
interactions these exponents need not to coincide with 
those of the RFIM. 
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